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The first year data from the Wilkinson Microwave Anisotropy Probe are used to place stringent 
constraints on the topology of the Universe. We search for pairs of circles on the sky with similar 
temperature patterns along each circle. We restrict the search to back-to-back circle pairs, and to nearly 
back-to-back circle pairs, as this covers the majority of the topologies that one might hope to detect in a 
nearly flat universe. We do not find any matched circles with radius greater than 25°. For a wide class of 
models, the nondetection rules out the possibility that we live in a universe with topology scale smaller 


than 24 Gpc. 
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Recently, the Wilkinson microwave anisotropy probe 
(WMAP) has produced a high resolution, low noise map 
of the temperature fluctuations in the cosmic microwave 
background (CMB) radiation [1]. Our goal is to place 
constraints on the topology of the Universe by searching 
for matched circles in this map. Here we report the results 
of a directed search for the most probable topologies. 
The results of the full search, and additional details about 
our methodology, will be the subject of a subsequent 
publication. 

The question we seek to address can be plainly stated: 
Is the Universe finite or infinite? Or, more precisely, what 
is the shape of space? Our technique for studying the 
shape, or topology, of the Universe is based on the simple 
observation that, if we live in a universe that is finite, light 
from a distant object will be able to reach us along more 
than one path. The one caveat is that the light must have 
sufficient time to reach us from multiple directions or, put 
another way, that the Universe is sufficiently small. The 
idea that space might be curled up in some complicated 
fashion has a long history. In 1900, Schwarzschild [2] 
considered the possibility that space may have nontrivial 
topology, and used the multiple image idea to place lower 
bounds on the size of the Universe. Recent progress is 
summarized in Ref. [3]. 

The results from the WMAP experiment [1] have deep- 
ened interest in the possibility of a finite universe. Several 
reported large scale anomalies are all potential signatures 
of a finite universe: the lack of large angle fluctuations 
[4], reported non-Gaussian features in the maps [5,6], and 
features in the power spectrum [7]. 

The local geometry of space constrains, but does not 
dictate, the topology of space. A host of astronomical 
observations supports the idea that space is locally homo- 
geneous and isotropic, so we may restrict our attention to 
the three-dimensional spaces of constant isotropic curva- 
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ture: Euclidean space E’, hyperbolic space H°, and 
spherical space $3. A useful way to view nontrivial to- 
pologies with these local geometries is to imagine the 
space being tiled by identical copies of a fundamental 
cell. For example, Euclidean space can be tiled by cubes, 
resulting in a three-torus topology. Assuming that light 
has sufficient time to cross the fundamental cell, an 
observer would see multiple copies of a single astronomi- 
cal object. To have the best chance of seeing “around the 
Universe” we should look for multiple images of the most 
distant reaches of the Universe. The last scattering surface 
(LSS), or decoupling surface, marks the edge of the 
visible Universe. Thus, looking for multiple imaging of 
the LSS is a powerful way to look for nontrivial topology. 

But how can we tell that the essentially random pattern 
of hot and cold spots on the LSS have been multiply 
imaged? The LSS is 2-sphere centered on the observer, 
and each copy of the observer will come with a copy of 
the LSS. If the copies are separated by a distance less than 
the diameter of the LSS, then the copies of the LSS will 
intersect. Since the intersection of two 2-spheres defines a 
circle, the two last scattering surfaces will intersect along 
circles. These circles are visible by both copies of the 
observer, but from opposite sides. Of course, the two 
copies are really one observer so, if space is sufficiently 
small, the cosmic microwave background radiation from 
the LSS will have patterns of hot and cold spots that 
match around circles [8]. The key assumption in this 
analysis is that the CMB fluctuations come primarily 
from the LSS and are due to density and potential terms 
at the LSS. 

Implementing our matched circle test is straightfor- 
ward but computationally intensive. The general search 
must explore a six-dimensional parameter space. The 
parameters are the location of the first circle center 
(01, $i), the location of the second circle center 
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(82, p2), the angular radius of the circle a, and the rela- 
tive phase of the two circles 6. We use the HEALPIX [9] 
scheme to define our search grid on the sky. A resolution r 
HEALPIX grid divides the sky into N = 12N2,, equal area 
pixels, where Ngige = 2”. Similar angular resolutions were 
used to step through @ and £. Thus, the total number of 
circles being compared scales as ~N*, and each compari- 
son takes ~N!/? operations. The simplest implementation 
of the search at the resolution of the WMAP data (r = 9) 
would take greater than 10% operations. This is not 
computationally feasible. However, with the algorithms 
outlined below, we can carry out a complete search for 
most likely topologies using the WMAP data. 

The WMAP data suggest that the Universe is very 
nearly spatially flat, with a density parameter Qg = 
1.02 + 0.02 [4]. Our universe is either Euclidean, or its 
radius of curvature is large compared to the radius of the 
LSS. For topology to be observable using our matched 
circle technique, we require that the distance to our 
nearest copy is less than the diameter of the LSS. This 
in turn implies that, near our location and in at least 
one direction, the fundamental cell is small compared 
to the radius of curvature. Given the observational con- 
straint on the curvature radius, there are strong con- 
straints on the types of spherical topologies that might 
be detected [10], and on where we can be located in a 
hyperbolic manifold for the topology to be detectable 
[11,12]. Naturally, the near flatness of the Universe places 
no restrictions on the observability of the Euclidean top- 
ologies. Remarkably, the largest matching circles in most 
of the topologies we might hope to detect will be back to 
back on the sky or nearly so. This immediately reduces 
the search space from six to four dimensions. This result 
holds for most circle pairs in nine of the ten Euclidean 
topologies (the exceptions are the largest circle pairs 
formed by twisted identifications, and the next-to-largest 
circle pairs in the Hantzsche-Wendt space). The result is 
less obvious in spherical space [11], but is nonetheless 
exact for a large class of such spaces (single action mani- 
folds). All others (double or linked action manifolds) 
predict slightly less than 180° separations between the 
circle centers for generic locations of the observer, but 
larger offsets for special locations. For topology to be 
detectable in hyperbolic manifolds, we would need to live 
near a “cusp” or “horn,” where the topology is of the 
form R X T? or R? X S [12] and the matching circles are 
back to back on the sky or nearly so. The possibility of a 
slight offset motivates the final search described in this 
Letter, where we consider circles separated by between 
170° and 180°. The possibility of a less probable topology 
requires the full search which will be described in a 
future publication. 

How do we compare circles? Around each pixel i, we 
draw a circle of radius «æ and linearly interpolate values at 
n = 2'*! points along the circle. We then Fourier trans- 
form each circle: T;(@) = X „Tim exp(im@) and compare 
circle pairs, with equal weight for all angular scales: 
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The i, j label the circle centers. With this definition S; a 
1 for a perfect match. For two random circles, the expec- 
tation value is zero, (S;;) = 0. We estimate (S?) below. 
We can speed the calculation of the S statistic by re- 
writing it as S;(a,B)= > mSme "8, where sm = 
2MT in TAX an Tal? + |T;,|?)]. By performing an in- 
verse fast Fourier transform ‘of the sm, we get S;;(@, B) at 
a cost of ~N'/? logN operations. This reduces the cost 
of the back-to-back search from ~N5⁄ to ~N? logN. 
We further speed the calculation by using a hierarchical 
approach: We first search at r = 7 and identify the 5000 
best circles at a fixed value of a. We then refine the 
search around these best pixel centers and radii and 
complete the search at r = 9. Smax(@) is the best match 
value of S;;(a@, £) for all pairs of circles within a specified 
range of angular separations, extremized over all relative 
phases £. 

To test the search algorithm, we generated a simulated 
CMB sky for a finite universe with cubic three-torus 
topology with side length 0.513 in units of the radius of 
the LSS. When only the Sachs-Wolfe effect [13] was 
included, the algorithm found nearly perfect circles: 
Smax = 0.99 with pixelization effects accounting for the 
residual errors. Real world effects significantly degrade 
the quality of the circles. The dominant sources of noise 
are velocities at the LSS. Following the approach outlined 
in the appendix of [14], we generate CMB skies with 
initial Gaussian random amplitudes for each k mode 
and then compute the transfer function for each value 
of k. This approach implicitly includes all of the key 
physics: finite thickness LSS, reionization, the integrated 
Sachs-Wolfe (ISW) effect, and the Doppler term. Since 
topology has little effect on the power spectrum for 
multipoles / > 20, we use the best fit parameters based 
on the analysis of the WMAP temperature and polariza- 
tion data [4]. The ratio of the amplitude of the potential 
and density terms (which contribute to the signal for the 
circle statistic) to the velocity and ISW terms (which 
contribute to the noise for the circle statistic) depends 
weakly on the basic cosmological parameters: Within the 
range consistent with the WMAP data, the variation is 
small. For the simulated sky, we assume that the primor- 
dial potential fluctuations had the usual 1/k* power spec- 
trum. The simulation included detector noise [15] and 
the effects of the WMAP beams [16]. The simulation 
did not include gravitational lensing of the CMB as the 
lensing deflection angle in the standard cosmology is 
small [17], ~6’, compared to the smoothing scale used 
in the analysis, 40’. When the search was performed on 
these realistic skies, the quality of circle matches was 
degraded. Figure 1 shows the results of the search on the 
simulated sky: The peaks in the plots correspond to radii 
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FIG. 1 (color online). The maximum value of the circle 
statistic as a function of radius, œ (in degrees), for a simulated 
finite universe model. The peaks in the plots correspond to 
positions of matched circles. The dashed (blue) line corre- 
sponds to the detection threshold described below. 


at which a matched circle was detected. The largest 
circles had Smas = 0.75, and the best match value 
decreased for smaller circles. The pure Doppler 
terms are correlated for large circles (@ > 45°), and 
anticorrelated for small circles. The code was able to 
find all 101 predicted circle pairs with radii greater than 
25°: The poorest detection had an Smax value of 0.45. 
While we only simulated a three-torus, the degradation of 
the circle match as a function of radius will be nearly 
identical for other topologies as it is determined by local 
physics at angular scales well below the circle radius. We 
will report in more detail on such tests in a future 
publication. 

For the circle search, we generated a low foreground 
CMB map from the WMAP data. Outside the WMAP Kp2 
cut, we used a noise-weighted combination of the WM AP 
Q, V, and W band maps. Using the template approach of 
[18], these maps were corrected for dust, foreground, and 
synchrotron emission. This map was smoothed to the 
resolution of the Q band map. Inside the Kp2 cut, we 
used the internal linear combination map [18]. 

The false positive signal can be estimated by approx- 
imating the CMB sky as a Gaussian random field 
with coherence angle 0.~ 0.7°. There are roughly 
Neire = 277/02 independent circles on the sky of radius 
a. Along each circle, there are 277 sin(a@)/0, independent 
patches and also 27rsin(a)/0. independent orienta- 
tions. Thus, the back-to-back search considered 
Necarch = 477 sin(a)/@2 possible circle pairs. If we 
treat the patches as independent, then (S?)!/? = 
1/./2asin(a)/6.. The number of circles expected 
above a given threshold Sọ is M(S >So) = 
Neearch/2erfe(So//2(S7)). Thus, inverting the comple- 
mentary error function yields the maximum value for 
any false positive circle at each radius: 
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This estimate is sensitive to 0., and excludes correlations 
between circles. 

We use ‘‘scrambled”’ versions of the WMAP cleaned 
sky map to obtain a more accurate estimate of the false 
positive rate. We generate scrambled versions of the true 
sky by taking the spherical harmonic transform of the 
map and then randomly exchanging a,,, values at fixed l. 
This scrambling generates new maps with the same two- 
point function but different phase correlations. The Snax 
value for the best fit back-to-back circle found at each 
radius is plotted as a function of «œ in Fig. 2. The false 
positive signal is well fit by a slightly modified form of 
the analytical estimate (2). Based on the analytical esti- 
mate, we set a detection threshold so that fewer than 1 in 
100 random skies should generate a false match. This 
threshold is shown in Figs. 1-3. 

Figure 3 shows the result of the back-to-back search 
on the WMAP data. The search checked circle pairs that 
were both matched in phase (points along both circle 
labeled in a clockwise direction) and flipped in phase. 
These two cases correspond to searches for nonorientable 
and orientable topologies. We did not detect any statisti- 
cally significant circle matches in either search. The 
simulations show that the minimum signal expected 
for a circle of radius greater than 25° is Smax = 0.45, 
which is above the false detection threshold. Since we 
did not detect this signal, we can exclude any universe 
that predicts circles larger than this critical size. In a 
broad class of topologies, this constrains the minimum 
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FIG. 2 (color online). The maximum value of the circle 
statistic as a function of radius for back-to-back circle pairs 
on scrambled CMB skies. The results are shown for seven 
different realizations. The white line is a fit for the false 
positive rate. The dashed (blue) line shows the false detection 
threshold set so that fewer than 1 in 100 simulations yield a 
false event. 
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FIG. 3 (color online). The maximum value of the circle 
statistic as a function of radius for the WMAP data for back- 
to-back circles. The dotted (red) line in the figure is for an 
orientable topology. The solid line is for nonoriented topolo- 
gies. The white line is the expected false detection level. The 
dashed (blue) line is the detection threshold. The spike at 90° is 
due to a trivial match between a circle centered at one point and 
its copy centered around the antipodal point. 


translation distance d = 2R,atan[tan(R;ss/ Re) COS(Qmin) |, 
where R. is the radius of curvature, and R;,, is the 
distance to the last scattering surface. Note that one 
recovers the Euclidean formula for R, — œ and the hy- 
perbolic formula for Re — iR,. For a flat or nearly flat 
universe, models with d < 2R),,cos(25°) = 24 Gpe are 
excluded. Note that the search excludes the Poincaré 
dodecahedron suggested in [7] as this model predicts 
back-to-back circles of radius 35°. 

By considering a broader set of circles, we can extend 
the search to other possible topologies. Figure 4 shows the 
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FIG. 4 (color online). The maximum value of the circle 
statistic as a function of radius for the WMAP data for circles 
whose centers are separated by greater than 170°. The solid 
line in the figure is for an orientable topology. The dotted (red) 
line is the expected false detection level and the dashed (blue) 
line is the detection threshold. 
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result of a search for nearly back-to-back circles: circles 
whose centers are separated by more than 170°. The 
restricted search requires searching far more circles 
than the back-to-back search, which increases the false 
positive signal. This increases the minimum detectable 
angular radius to 28°. The signal at a ~ 85°—90° is due to 
nearly overlapping and osculating circles. This feature is 
also seen in scrambled skies. 

Has this search ruled out the possibility that we live in 
a finite universe? No, it has only ruled out a broad class of 
finite universe models smaller than a characteristic size. 
By extending the search to all possible orientations, we 
will be able to exclude all topologies out to ~24 Gpc. 
More directed searches can extend the result for specific 
manifolds somewhat beyond 24 Gpc. With lower noise 
and higher resolution CMB maps (from the extended 
mission of WMAP and from Planck), we will be able to 
search for smaller circles and extend the limit to 
~28 Gpc. If the universe is larger than this, the circle 
statistic will not be able to constrain its shape. 
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